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Abstract 

For closed-shell systems, the local density approximation (LDA) and the LYP, BLYP, and B3LYP 
functionals are shown to be compatible with reference-state one-particle density-matrix theory, 
where this recently introduced formalism is based on Brueckner-orbital theory and an energy 
functional that includes exact exchange and a non-universal correlation-energy functional. The 
method is demonstrated to reduce to a density functional theory when the exchange-correlation 
energy-functional has a simplified form, i.e., its integrand contains only the coordinates of two 
electron, say ri and r2, and it has a Dirac delta function 6{ti — Y2) as a factor. Since Brueckner 
and Hartree-Fock orbitals are often very similar, any local exchange functional that works well 
with Hartree-Fock theory is a reasonable approximation with reference-state one-particle density- 
matrix theory. The LDA approximation is also a reasonable approximation. However, the Colle- 
Salvetti correlation-energy functional, and the LYP variant, are not ideal for the method, since 
these are universal functionals. Nevertheless, they appear to provide reasonable approximations. 
The B3LYP functional is derived using a linear combination of two functionals: One is the BLYP 
functional; the other uses exact exchange and a correlation-energy functional from the LDA. 

PACS numbers: 



I. INTRODUCTION 



The Kohn-Sham version of density functional theory (DFT)_plays a major role in both 



quantum chemistry and condensed matter physics 
approximation (LDA) has been widely used 



sory IJJI* i ipiays a 



most successful functional, a hybrid one 



S 
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7|. The local density 
br the solid state, while for molecules, the 
is known as B3LYP js, 13|, where this 



functional contains three parameters, two correlation-energy functionals, the Dirac-exchange 
: unctional with the a correction, and exact exchange. The LYP correlation-energy functional 
lJ | - a kevcomponent of B3LYP - is derived from the CoUe-Salvetti correlation-energy 
functional [15[- 

Recently, a density-matrix formalism was developed called reference-state one-particle 
density-matrix theory 16|, |l7|. This method employs generalized Hartree-Fock equations 



containing the exact exchange-potential and a nonlocal correlation-potential, where these 
equations are obtained using the Brillouin-Brueckner condition (for non-variational energy 
functionals) or functional minimization (for variational functionals). The method uses en- 
ergy functionals that yield the target-state energy when the reference state - or its cor- 
responding one-particle density-matrix - is constructed from Brueckner orbitals. Unlike 
traditional DFT approaches, the ti-representable problem does not appear in the method, 
nor the need to introduce functionals defined by a constrained search. The correlation- 
energy functionals in the method are non-universal, in the sense that they depend on the 
external potential. Nevertheless, model systems can still be used to derive universal energy- 
functionals. For example, the Colle-Salvetti functional was shown to be compatible with 
the method as well as a one-particle-density-matrix variant of the LDA. In addition, using 
time-independent many-body perturbation theory 19 1, diagrammatic expansions were 



given for the non-variational correlation-energy functionals that are expressed in terms of 
orbitals and orbital energies. When restrictions are placed on the orbital energies, the in- 
dividual diagrams are shown to explicitly depend on the one-particle density-matrix of the 
reference state (and the external potential). 

Below we demonstrate that for systems that can be treated with a simplified exchange- 
correlation energy-functional - leading to a local exchange-correlation potential - the 
reference-state one-particle density-matrix formalism reduces to a density functional for- 
malism, where the density that appears within the method is the one from the reference 



2 



state, and not the target state. This approach differs in this way and in others from the 
formahsm by Lindgren and Salomonson 2^, where their method is a generahzation of the 
Hartree-Fock Kohn-Sham approach 21| and yields orbitals that are believed to be very 
similar to Brueckner orbitals. 

The DFT functionals derived from a uniform electron gas are also shown to be compatible 
with the reference-state one-particle density-matrix theory, where this compatibility appears 
because the density of the target state for a uniform electron gas is the same as the density 
of the reference state and, also, the plane-wave states are Brueckner orbitals. In additional, 
functionals that approximates the exact exchange-energy can also be used, since the exact 
exchange-functional for Brueckner orbital theory, Hartree-Fock theory, and DFT are the 
same. Furthermore, functionals derived using the one-particle density-matrix of the reference 
state are compatible, assuming the errors due to the external potential are small - if the 
functional is a universal one. Using these observations, the LDA, and the BLYP and B3LYP 
functionals are shown to be compatible with the reference-state one-particle density-matrix 
formalism. This puts these functionals within an alternative DFT approach that avoids, for 
the most part, existence theorems and constraint searches. 

II. REFERENCE STATE ONE PARTICLE DENSITY MATRIX THEORY 

We seek an eigenfunction |\E^) on the time- independent Hamiltonian operator: 

ij i jkl 



where the integrals are written using chemist's notation 

Mkl] = Y.j ^;(xi)^,(xi)rr2Vfc(x2)V'Kx2) dv.dv,, (2) 

and the spatial and spin coordinates, r and oj, are denoted collectively by x. The eigenfunc- 
tion of interest, or target state |\E'), yields the electronic energy £: 

Consider a single-determinantal reference-state, say 1$), where this state is a first-order 
approximation of the target state; the first-order energy is 

= ($|//|$). (4) 



Adding and subtracting on the left side of Eq. Q, gives 

8 = Ei + {8-E^) = E^[^ + ^eo[$], (5) 

where this Eq. defines the correlation energy £^co['^'] as {8 — Ei). 

There is a one-to-one correspondence between determinant states and their one- par ticle 
density-matrices 0,12^, where these density- matrices 7 are given by 0, 2^, 26, 27, 3| 



7(x,x')= Yl ^-WC(x'); 7"l<^>), (6) 
weji/'o— >*} 

and the notation on the right-hand-side indicates the one-to-one correspondence; further- 
more, the sum over w includes only the occupied orbitals of the reference state |$): This 
orbital set is denoted by {i/Jo the unoccupied set is denoted by {ipu Un- 

less stated otherwise, two sets of occupied or unoccupied orbitals that differ by a unitary 
transformation are considered equivalent. 

Because of the one-to-one correspondence mentioned above, any functional that depends 
on the reference state |$) can be transformed into one that depends on the one-particle 
density-matrix 7. For example, by requiring the external potential v{r) to be a spin- free 
operator, the first-order energy can be written as 

i?i[7] = J [-iVh(x,x')]^,^^rfx + J t;(r)7(x,x)dx + Ej[7] + i?j7], (7) 

where the Coulomb and exchange energies have their usual forms: 

Ejb] = \JJ ^r2S(xi,xi)7(x2,X2)dxi(ix2, (8) 

-E^l] = \jj ^r2S(xiiX2)7(x2,xi)rfxi(ix2, (9) 

and, henceforth, an integration over x implies an integration over r and a summation over u. 

A. Hartree Fock Theory 

Let the Fock operator be denoted F^, the functional derivative of the first-order energy 
Eil'y] gives the kernel of this operator 

F,i^u^2) = t4^^ (10) 
57(x2,xi) 

= (5(x2-xi) l^-iVa +t;(r2) + j r23S(x3,X3) dxg^ + t;;^ (xi , Xg) , 



where the two-body function, f^(xi,X2), is obtained from the exchange energy: 

^^x(Xl,X2) = ^- = -ri2 7(xi,X2), (11) 

d7(X2,Xi) 

and this function is the kernel of the exchange operator, denoted by v^. 

Minimizing the functional i?i[7], subject to the constraint that the one-particle density- 
matrix comes from a single-determinantal state, yields -E'i[t], where f is the one-particle 
density-matrix of the Hartree-Fock reference-state |G): 

f(x,x')= ^-WC(x'), (12) 

and this density-matrix satisfies the following condition: 

{ipr\Ff\'ipJi = ^Ijw ^ bl^o ^ r], tpr ^bPu^r}, (13a) 

and the right-hand-side notation indicates that -i/^^ is a member of the set of occupied orbital 
that determines f (and |0)); 'ipr is unoccupied. The operator form of this Eq. is 

(i - r) Frf = 0, (13b) 

where the identity operator 1, using a Hartree-Fock basis set, is given by 

i = Yl + E \^r){i^r\. (14) 

A unique set of occupied and unoccupied orbitals is obtained by requiring the occupied 
and unoccupied blocks of Ff to be diagonal: 

F,rM = eirM, tfjl e {iJo f}, (15a) 

F^^;(x) = £;v^;(x), e {^u ^ t}. (15b) 

The orbitals sets that satisfy Eqs. (fT^ and (fT3j) are denoted by {tpl ^ ^, -^f } and 
{ipl ^ T, F-f }, indicating that they are uniquely determined by f and Ff. 
Substituting Eq. (fTUI) into Eqs. (fT^ gives the canonical Hartree-Fock Eqs: 

(^-iV? + i;(ri) + j rf^^f (X2,X2) rfx2 + v'M^ Va'(xi) = ^[^[(xi). (16) 
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B. Variational Brueckner orbital theory 



Brueckner orbital theory 



29 



30 



31 



32 



33, 



35 



39 



40 



41 



42 



Q is 



a generalization of Hartree-Fock theory that utilizes a single-determinantal state that has 



the maximum overlap with the target state 



45l | . We now review a variant of Brueckner 



orbital theory that is used to derive the reference-state one-particle density-matrix theory 



In this previous work, we introduced four trial wavefunctions - say 1^1/$'), where {rj = 
I, II, III, and IV) - that are defined with respect to a single-determinantal reference-state 
1$). The first trial-wavefunction is simply the target state of interest, with the 

single excitations removed: 

l^g') = (i-At)l^), (17) 

where the projector for the singly-excited states is 



and the singly-excited state are 



(18) 



(19) 



The subspace is completely determine d by 1$); is also invariant to a unitary trans- 
formation of occupied, or virtual, orbitals 36|. 

The second trial-wavefunction is defined with respect to the target state expressed 

by an exponential ansatz: = e'^*|<l')), where is generated by removing the single- 

excitation amplitudes Sf from the cluster-operator S: 



1$). 



(20) 



The third trial-wavefunction ') can be generated by its wave-operator: 



(21) 



that can be expressed in an exponential form: {Clq, = 6"^*!$)), where can be written as 
a sum n-body excitations, with the exclusion of a one-body operator: 



(22) 



The wave operator IS ClS db solution to the coupled cluster equations 
with the single excitation portion removed: 



49 



50 



51 



18|, 



IS 



(1 - p*) (hq^) = 0, 

V / op,cn 



(23) 



where only the open (op) and connected (cn) portions enter into the relation. This expression 
defines the trial functional |\E'^"') using Eq. pij) . 

The fourth trial wavefunctions I^E'^^') is not considered here, exce pt t o mention that it is 
obtained by solving the configuration-interaction equations jloL l22l . l54| in an approximate 
way, i.e., by neglecting the single-excitation portion. 

All of the trial states 1^1'$'') share the property that they contain no single excitations, i.e., 
(P*|\l'^') = 0), and they generate the target state l^') when their reference state satisfies 
(1$) = |0)), where |G) is the determinantal state constructed from occupied Bruckner 
orbitals. In other words, we have 

1^0 ) = 1^)- (24) 
Because of the one-to-one correspondence mentioned above, the trial wavefunctions de- 
pend on the one-particle density-matrix; so we can write |\E'iy'''). Using this notation, we can 
define variational energy-functionals that depend on the one-particle density-matrix: 



(25) 



where the last relation defines the correlation-energy functional i?co [7] as (-E,y[7] — -E-iH)- 
Substituting Eq. ((Tj) into Eqs. (0) and gives the following: 



^ = / [-lVh(x,x')]^,^^rfx + I t;(r)7(x,x)dx + Ej[7]+^.j7], 



(26a) 



I [-iV^7(x,x')]^,^^dx + I t;(r)7(x,x)dx + i5;j[7] + ^^,'[7], (26b) 



where the exchange-correlation energy and exchange-correlation energy-functionals are, re- 
spectively, defined by 



^xc[7] = ^x[7]+^co[7], 



(27a) 
(27b) 



We now generalize many of the Hartree-Fock relations from the previous subsection 
simply by replacing i?i[7] with -^^[7] and the Fock operator with generalized, or exact, 
Fock operators, C*'''. 



The functional derivative of our energy functionals yield two-body functions that serve 
as the kernels of the exact Fock operators: 

Cr(xi,X2) = (28) 

= - xi) (^-iV^ + v{r2) + j r2t7(x3,X3) d^^ + (xi, X2), 

where the kernels of the exchange-correlation operators, t'^^'(xi, X2), are obtained from the 
exchange-correlation energy-functionals: 

^x\^(xi,X2) = = - rr2S(xi,X2), (29) 

57(x2,xi) 57(x2,xi) 

where the last relation uses Eqs. ()27b|) and (fTTj). 

We now require the target state to be a ground state; using the variation theorem and 
Eq. ()24|). the minimizing of the functionals -^^[7] - subject to the constraint that the one- 
particle density-matrix comes from a single-determinantal state - yields 

S = E,[Ti (30) 

and from Eqs. (j^UI) and (fTTj) we have 

^eo[r] = (31) 

^xcM = mr], (32) 

where r is the one-particle density- matrix of the Brueckner reference-state |6): 

r(x,x')= V^-(x)C(x'), (33) 

and the Brueckner orbitals satisfy the following equivalent conditions: 

{'iIjXC;'Mu,) =0; {iIjo ^ r}, 7/;,, e ^ r}, (34a) 

(i - r) eV = 0, (34b) 

where these orbitals do not depend of - any trial wavefunction gives the same results. 

A unique set of occupied and unoccupied orbitals is obtained by requiring the occupied 
and unoccupied blocks of C,^^ to be diagonal: 

C|''V:(x) = CC(x), C e {^o ^ r}, (35a) 
CV;(x) = C^;(x), e {^n ^ r}. (35b) 



Henceforth, the orbitals sets that satisfy Eqs. (jH^ and (jSHl) are denoted by {ip'^ ^ t,Ct ^} 
and {ipl^ ^ T, Cr''}, indicating that they are determined by r and CI'''- Since theses orbitals, 
and their energies, can, perhaps, depend on t], it is more precise to denote then by ipj^ and 
^j^**, but we suppress the t] superscripts to keep the notation less cluttered. 

Substituting Eq. into Eqs. (jH^ gives generalized, canonical Hartree-Fock Eqs: 



iVl + v{ri) + j rr2V(x2,X2)dx2 + z>-(xi)^ V^^xi) = CV'Kxi). 



(36) 



III. CLOSED-SHELL CASE 

Our interest here is in cases where the Hamiltonian is spin-free - it contains no spin 
coordinates - and the target state is a ground-state singlet that is well described by a 
closed-shell reference-state |$). In these cases - and many others - it is appropriate to use 
a set of spatially restricted orbitals, given by 

i'jcri^) = Xi{.r)a{uj)] a = a,p, (37) 

so that 1$) is determined by a set of doubly-occupied spatial orbitals, denoted by {xo^ 
where this set also determines the virtual set, denoted by {xw^ Two sets of orbitals 
that differ by a unitary transformation are, again, considered equivalent. 
By definition, the spinless, one-particle density-matrix is given by 3, |3| 

Pi(ri,r2) = ^7(ri,c^;r2,w), (38) 

and from Eqs. pTjl and © with {w = xa) we have 

pi(ri,r2)= 2^ Xx(ri)x:(r2); Pi " |<f). (39) 

By using the restricted orbitals defined by Eq. (jSZj), it is easily demonstrated that the one- 
particle density-matrix 7 is determined by the spinless one, as indicated by the following 
relation: 

7(xi,X2) = ^Pi(ri,r2)5a;i<^2- (40) 

Hence, any functional of 7 now becomes a functional of pi. In particular, for a spin- free 
Hamiltonian, the spin- variable summations are easily performed, yielding explicit functionals 
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of Pi. For example, the electronic energy and our energy- functional, Eqs. © and ()25p. can 
be written as 

£ = Ei[pi]+S,o[pi], (41a) 

E,[p,] = E,[p,] + EM, (41b) 

or the representations of these functionals, given by Eqs. (j^ . reduce to 

^ = J [-^V^Pi(r,r')]^,^> + J vir)pir)dr + Ej[p]+£^,[pi], (42a) 

Er,[Pi] = J H'^l Pi(r, r')]^,^dr + J t;(r)p(r) dr + + E^^[p,], (42b) 

where the definitions given by Eqs. (|77|) can be expressed using a dependence upon pi. 

£^c[pi] = E^Ipi] + £co[pi], (43a) 
E^J[p^]=E^[p,] + EM■, (43b) 
furthermore, the first-order energy, Eq. ((7j), becomes 

Ei[pi] = J [-iV^pi(r,r')]^,^/r + J viv)piv) dv + Ej[p] + E^Pi]. (44) 
where the Coulomb and exchange energies, Eq. (jH)) and (P), reduce to 

Ej[p] = \jj ri2p{ri)p{r2)dridr2, (45) 
-E^IPi] = \jj r^ipiiri,T2)piir2,ri) dridr2, (46) 
and the particle density is 

p(r) = pi(r,r). (47) 

A. Hartree Fock Theory 

Using Eqs. (jllj) and ()4U|) it is readily observed that the kernel of the exchange operator 
is given by 

z;;j(xi,X2) = <i(ri,r2)(5<^,,^2, (48) 

where 

<^(ri,r2)= -irr2Vi(ri,r2). (49) 
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Therefore, for an arbitrary function, say 0, we have 

i);^</>(xi) ={)^i0(xi), (50) 
where -Uxl^ii^a) and t;^i(ri,r2) serve as the kernels of v1 and v'^, respectively: 

^l(t){^l) = j t'x(xi,X2)0(x2) rfX2, (51) 

{)^V'(ri) = y"f^i(ri,r2)0'(r2) dr2. (52) 
Using Eqs. (pSjl and (jlUI) . the kernel of the Fock operators, Eq. ()10|) . can be written as 

F^(xi,X2) = 5^i<^2Fpi(ri,r2), (53a) 

where 

Fp,(ri,r2) = 5(r2-ri) (^-1^^^ + v{v2) + j r:,^ p{v:i) dv^ +<Kri,r2), (53b) 

Substituting Eqs. (|37j) and ()53a|l into (fT5|) . and summing over ciJ2 we get 

y F^,(ri,r2)xf (r2) dv^aiuo^) = 5f xf (ri)a(^i). (54) 

By letting F^^(ri, r2) serve as the kernel of the spin-free Fock-operator F^^, we have 

4xf (ri) = efxf (ri), (55) 

where is the spinless one-particle density-matrix of the Hartree-Fock state: 

f(xi,X2) = ^Pi(ri,r2)5^ia;2> (56) 

and the Hartree-Fock spatial-orbitals are denoted in an analogous way as the spin-orbitals 
from Sec. pi A|) . e.g., replaces f. 

Substituting Eqs. ()53b|) into and using gives 

(-iV? + v{v{) + j t^^q{v^) dr, - {)|^(ri)) xf (ri) = efxf (ri). (57) 

Eqs. ()55|1 and ()57|1 are the closed-shell spin-free forms of Eqs. (fTH|) and (fTI)|) : the spin- free 
forms or Eqs. ()10p and are easily proven to satisfy the following relations: 



where, as mentioned previously, these functions are given by Eqs, ()53b|) and (p9|). respec- 
tively. 
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B. Brueckner orbital theory 



We now generalize many of the Hartree-Fock relations from the previous subsection 
simply by replacing the exchange operator f)^ with the exchange- correlation operator z>^^. 
For the closed shell case under consideration, it can be easily demonstrated that the 

n 

kernels of the exchange-correlation operators are given by [55|| 

^x7(xi,X2) = C(ri,r2)5^„.,. (60) 

Therefore, for an arbitrary function, say 0, we have 

/>-0(xi) = j>^f0(xi), (61) 

where z/^^(xi,X2) and zy^J''(ri, r2) serve as the kernels of i>^^ and P^^, respectively. 

Using Eqs. and pn|) . the kernel of the exact Fock operators, Eq. (|^. can be written 

as 

C;'"(xi,X2) = 5^,^,C^^'(ri,r2), (62a) 

where 

C;?(ri,r2)=5(r2-ri)(^-iV,2 + t;(r2)+ J r^t^rg) (irg^ + C(ri, r2). (62b) 
Substituting Eqs. (jTTj) and ()62a|l into and summing over uj2 we get 

I C<f(ri,r2)xf (r2) dr.aM = ifxf{r,)a{uj,). (63) 
By letting C^'^'(ri,r2) serve as the kernels of the exact, spin-free Fock-operator we have 

C'xf (ri) = efxf (ri), (64) 
where Qi is the spinless one-particle density of the Brueckner state: 

r(xi,X2) = ^£)i(ri,r2)5c.ia;2, (65) 
and the Brueckner spatial-orbitals are denoted in an analogous way as the spin-orbitals from 



Sec. piB|l . e.g., qi replaces r. 

Substituting Eqs. (l62H|l into (jHSl) gives 



(^-\Vl + v{r,) + j r~,^Q{r,) dv, + i)-''(ri)^ xf (ri) = ^xf (ri). (66) 
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Eqs. and (jHEl) are the closed-shell spin-free forms of Eqs. and (j^: the spin-free 
forms or Eqs. pH)) and (I^Uj) can be shown to satisfy the following relations j55|: 

- '''' 

.™(n,r.) ^ iM^. (68) 

dpi(r2, ri) 

Substituting Eq. ()43bj) into and using (jKHj) . we have 

z.^J''(ri,r2) = <Kri,r2) + /^eT(ri,r2), (69) 

where the kernel of the exchange potential v'^'^{vi, T2) is given by Eq. ()49|): the kernel of the 
correlation potential u^^{yi,Y2) is defined by 

The operator form of Eq. ()69p is given by 

= + ^cT- (71) 



IV. APPROXIMATIONS 



The target state |\E') can be generated from a reference state |$) by its wave-operator fi^: 

1^) = 17^1$) = |^^) + |^^e), (72) 

where we have partitioned the target state into two terms: |\E'$) and |\E'<j, 0). Consider 
the diagrammatic perturbation-expansion for |\E'<i),e) is then defined as the sum of 

all diagrammatic terms from the perturbation expansion that must vanish for any external 
potential when the reference state is the Brueckner one: 

|*e,e) = 0. (73) 

Hence, 1^^$) is a trial wavefunction that fulfills all of the requirements of the ones considered 
in Sec. Ill Bl furthermore, by our definition, 1^^$) is the best trial wavefunction, since it 
contains no extra terms that must vanish when (|$) = 

Consider a system with an external potential, say v', where its correlation energy Sco 
is known when the reference state is the Brueckner one with a corresponding one-particle 
density- matrix of r'; using Eq. (j!?T|) we have 

£co[r',v'] = E,4t',v'], (74) 
13 



where the correlation-energy functional from is denoted by Eco and we have indicated 
that these functionals depend on the external potential. Hence, a reasonable approximation 
for Eco at other potentials v and one-particle density-matrices 7 is given by 

Ecob,v] = Sco[T',v\^,=^y=^y (75) 

Now consider the four trial wavefunctions which can also be obtained using a wave oper- 
ator: 

\¥^') = n'^'\<i>), (76) 

where from the previous notation we have (fi^"' = Q^). Consider the first two cases: 
(?7 = I and II) , where the wave-operators expansion are easily obtained by considering Lind- 
gren's formalism By examining this expansion, it is readily seen that the trial 



wavefunctions contain terms, say I^I/^'e), that, by necessity, vanish when (|$) = |6)): 

Hence, even though the approximation (fTSj) is reasonable for [r] = 1 and II): 

Kl'b,v] = Sco[r',v']^r'=,y=v), (78) 

it neglects the contribution to Ej^^'' arising from 

The wave operator diagrams for [r] = III) are more complicated to describe and are given 
elsewhere jsSj. However, it is obvious from the previous work \y\ that the term I^I/^^e) is 
considerably smaller in magnitude than for (r/ = I and II) , suggesting that the approximation 
flTSj) is reasonable, where this type of approximation was considered previously [l^. It is 
quite possible that l^l^^e) is zero, so, in that case (|\1/^"') = I^P*)), this was previously 
presumed to be true |l7|, but it was not discussed in detail. 

The wave-operator diagrams for [r] = IV) contain unlinked terms and will not be consid- 
ered, except to mention that I^&^^q) is nonzero. 

A. Local Density Approximation 

The correlation energy of an electron gas - denoted by ^^r"^ ^ is given by the correlation 
energy Sco for a system where the external potential v is constant. Since the orbitals for an 
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electron e; Brueckner, and using Eq. ()31|) . we have 

£ir'{rg] = £co[rg, v] = Eil^lTg, v], v = const, (79) 

where Tg is the one-particle density-matrix that is from the Brueckner reference-state of an 
electron gas, say |6g), where the Brueckner orbitals, defining \Qg), are given by plane waves 
when periodic boundary conditions are imposed; furthermore, the correlation energy Sco is 
unchanged if v is changed by a constant amount; i^cr°''[Tg] does not depend on v. (A constant 
shift in V does, however, shift the first-order energy Ei.) 

The electron-gas approximation for the third correlation-energy functional utilizes this 
functional, and it is given bv jl7|: 

En] ^ ^^r%]ir,=,), (80) 

and this approximation is a special case of Eq. (f7S|) . where there are no terms involving 
an external potential. (We can also use this approximation for E^o defined in the previous 
subsection.) 

For spatially restricted orbitals, given by Eq. ()37|) . the above Eq. becomes 

^r[pi] ^ ^r^b.](...=pO' (81) 

where gig is the spin-less, one-particle density-matrix of the Brueckner reference-state for 
an electron gas: 

rg(xi,X2) = ^^ig(ri,r2)(5^ic^2- (82) 

For a uniform electron gas, Eq. (|H^ is denoted by 

1 

7"ug(xi,X2) = -^'iug(ri,r2) 5^1,^2= (83) 

and the density that is from the Brueckner, one-particle reference-state of a uniform electron 
gas, say IG^g), is given by 

^ug(r) = ^iug(r,r), (84) 
and this density is constant, i.e., it does not depend on r; furthermore, we have 

^'ug = ^ug, (85) 
where n^g is the density from the target state, say l^&ug), of a uniform electron gas: 

l^^l) = l^ug), (86) 
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where we have used Eq. 

The correlation energy S^'^^^ of a uniform electron gas depends on the constant density 
n-ag] it is given by the following relation: 

£:rKg) = iVe-Kg), (87) 

where eJJ^ is the correlation-energy per particle and N is the number of particles. Using 
Eq. (jH3j) gives 

^r(f?ug) = A^C(^ug). (88) 

Following Kohn and Sham Q], but using Eq. (jHHj) instead of (|H7j). we construct the following 
functional for an electron gas: 

^'cT^g]^ j 69{r)eZ{gg)dr. (89) 

This functional can be expressed as one that appears to depend upon the one-particle 
density-matrix: 

^cr'K] - J J ^i9(r2,ri)5(ri - r2)e- (f?i,(r2, ri)) dv^ dv2, (90) 
where {gig{r,r) = gg{r)). Substituting the above expression into Eq. (jHTj) gives 

^rbi] ^ f f Pi(r2,ri)5(ri -r2)e-(pi(r2,ri)) dv, dv,. (91) 



The correlation potential z/^^'' is defined by its kernel, given by Eq. (jZOj). Because of the 
(5(ri— r2) factor, z>^i"' is a local operator when using the above expressions for Ej,™^; explicitly, 
using the results from Appendix^ we have the following approximate expressions: 

[Pi] ^ E:^[p] = I P(r)e-(p) dv, (92) 



where, for future reference, these approximations are denoted by E^^ and i/^^. 

Using the exchange-energy formula by Dirac for a uniform electron gas, the exchange 
portion can be treated in an identical manner as the correlation; this approach yields the 
following local exchange-operator and exchange energy: 

EM « E"/[p] = 5 (1) ' j p(r)-'/3 ar, (94) 
^'^-'l^l ^^pwV". (95) 



16 



where E^^ denotes the Dirac exchange-energy functional and v^^ is the corresponding ex- 
change potential. The above four expressions are identical to the ones for the LDA of DFT, 
except that - in the current approach - p is not required to be the density of an exact 
wavefunction |\E'); p is only required to be the density of a reference state |$), and this is 
assured to be the case. 



B. The Becke exchange correction 

For closed-shell systems that are well described by a single determinant, it is well known 
that the Bruckner and Hartree-Fock reference-states, |G) and |G), are usually very simi- 
lar. Hence, in these cases, the spinless, one-particle density-matrices from the Brueckner 
reference-state, gi, is similar to the one from the Hartree-Fock reference state, ^i. There- 
fore, an approximate exchange-functional i?x[pi] that depends on the density of the reference 
state, i.e., 

^xN~^x[p], (96) 

and is an accurate approximation with the Hartree-Fock approach, is also a reasonable 
approximation for the current approach. Of particular interest is the exchange functional 
obtained by Becke [57] - where this functional is an approximation of the Hartree-Fock 
exchange with empirical parameters fit to atomic systems; explicitly, this functional is given 

by 

EM = mp] + ^rip], (97) 

where AE^^^ is the Becke exchange correction to the uniform electron-gas term, E^^. Since 
this correction term contains a gradient dependence, the general form of the functional can 
be written using the following integral relation: 

EAp] = J Gx(p,Vp)dr. (98) 
By using Eqs. and and the following identity: 

Vip(ri) = [(Vi + V2)pi(r2, ri)]r,=ri ' (99) 
we obtain the following approximation: 

EMi]-K''[p] + ^r[p] = j I S{T2-r,)G, (pi(r2,ri),(Vi + V2)pi(r2,ri)) dndr2. 
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Using this expression for the exchange energy, the exchange potential is local; using the 
result from Appendix El and the definition for v^^ given by Eq. we have 

^ ^ ' ^ ^ ^ 5p(r) ^ V(r) ' ^ ^ 

where this expression defines v^^'^lr). 



C. The LYP Functional 

Consider the approximation given by Eq. ()78|) , where {rj = III) , and the correlation energy 
Sco - which is also a functional - is taken to be the one from the helium atom: 

where is the Brueckner, one-particle, density-matrix for the helium atom, and v^^ is the 
external potential for this system. Furthermore, since the Brueckner density-matrix r„c is 
approximately equal to the Hartree-Fock one, say f^^, we can write 

(102) 

The dominant portion of the correlation energy comes from the electron-repulsion contri- 
bution. The portion arising from the external potential is smaller; since v is a one-body 
operator, this term is probably treated well as part of the first-order energy Ei. Therefore, 
neglecting to make the substitution for the potential, i.e., {v^^ = v), should only yield a 
small error: 

~^co[r„.,^;J(f„,=^). (103) 
A well known approximation for Srn \t„ , - that is valid for other systems also - is given 



by Colle and Salvetti functional 
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^Co[r„e,^^„J ~^Co[^He]. (104) 

The above two relations suggest the following approximation: 

While this is probably a reasonable approximation, we do have some reservations about it. 
In particular, a universal functional is used to approximate one that depends on the external 
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potential. This violates, at least, the spirit of the method. We would therefore prefer to use 
the approximation given by Eq. ()101|) or ()102j) . Unfortunately, however, no such functionals 
are available. Therefore, let us proceed. 

For a closed-shell system, Eq. (jl(J5|) becomes 

^^o"'[Pi] ^ EM = ^:o1^iH.](,,,.=p,), (106) 

where ^iho is the Hartree-Fock spin-less one-particle density-matrix for the helium atom, and, 
for future reference, this approximations is denoted by E^l[pi]. An expression for E^l[pi] is 
derived by Lee, Parr, and Yang [l^ that is equivalent to the one derived by CoUe-Salvetti; 
it has the following form: 

ETo[pi] = j G'co(pi(ri,ri),Vi/9i(ri,ri),V?pi(ri,ri),(ViV2/9i(r2,ri))^^^^J dn, (107) 

and using this Eq. and the previous one, we have 

ET[pi] - j j (^(rz-rOGeo (pi(ri,ri),Vi/9i(ri,ri),V?pi(ri,ri),ViV2Pi(r2,ri)) dr^dv^. 

(108) 

When the functional derivative is taken to determine the kernel of the correlation potential 
z^^^"'(ri, r2) - as defined by Eq. (fTOj) - and £'co"^ is approximated by the above relation, the 
5{y2 — ri) term survives; the correlation potential P^^"' is a local operator. Denoting this 
approximation by z/^j'(r), we have 

^rw^^cTw- (109) 

Using the exchange potential given by Eq. ()100p - or any other local exchange potential 
- and the above local correlation operator, the exchange-correlation potential i)'^ - given 
by Eq. (|71|) - is also local. Therefore, the reference-state |$) obtained from the following 
non-interacting Hamiltonian: 

i i 

is also obtained from a local potential, where 

/vf (n) = v(v,) + / r^^Q^iv,) rfr2 + ^;-(r) + ^-(ri) + iP,!{v,), (111) 
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and Eq. (j66|) . for the occupied orbitals, becomes 

(-iV? + /vf (ri))x^nri) =Cx^nri). (112) 

If these equations are solved in a self consistent (SCF) manner, a series of spin-less, one- 
particle density-matrices are generated: p\, p\, pi, ■ ■ ■ ; for the nth iteration, we have 

(_iV? + ^;?(rO) xe\ri) = Xr(ri), (113) 

and this iteration yields a non-interacting state, say l^""*"^), with a one-particle density- 
matrix of Pi~^^; in addition, 1$"+^) is an eigenfunction of a non-interaction Hamiltonian, 

58|, each non- interacting states generated in 



Hs^ . From the Hohenberg and Kohn theorem 
the SCF procedure are unique functions of their electron densities, e.g., |$"(p)); furthermore, 
because of the one-to-one correspondence between one-particle density-matrices and non- 
interacting states, all one-particle density-matrices generated from the SCF procedure are 
determined by their density, i.e., p\{p)- And since for an SCF computation, these are the 
only one-particle density-matrices that are required, we can write 

EZipi] = EZ[p]. (114) 

Because of the general form of -E'^o[pi] given by Eq. ()107j) we need only an expression for 
(ViV2Pi(r2, ri))r2=ri- Following Lee, Parr, and Yang, we use the following approximation: 

(v,V2Pi(r2, rO) - \{?>n^r^p{v,f'' + ^ ^^'ff + ^V?p(rO, (115) 
V /r2=ri 5 36 p(ri) 4 

and we have 

EZ[p]^E'M (116) 

where i?'^P[p] is the LYP functional obtained by substituting the approximation given by 
Eq. (jllSp into the i?co[Pi] functional. Using the general of form of E^l[pi] given by Eq. ()107|) . 
we have 

E'ZnP] = J Ceo (p(r), Vp(r), VV(r)) dr. (117) 

Using integration by parts, Miehlich, Savin, StoU, and Preuss, eliminate the Laplacian terms 
and obtain the general form given by [s^ 

ElnP] = J Geo (p(r), Vp(r)) dr. (118) 
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Using Eqs. ()106|) . pi4|) . pi6|) . ()118p . and we obtain the following approximation for 
the correlation-energy functional: 

^ E'jJ'[p]= 1 1 5(r2-ri)ao(pi(r2,ri),(Vi + V2)pi(r2,rO) dr.dr^, 

(119) 

and using the results from Appendix the local correlation potential is given by 

^cr(r) ^ C^(r) = ^1^. (120) 

6p[r) 

Using these two approximations for E]^J'[p] and z>^i"^(r) yields an identical approach as the 
LYP method. If, in addition, the Becke exchange functional, given by Eq. pOO|) . is used, 
the method is equivalent to the one known as BLYP. 



D. The B3LYP Functional 

Using Eq. ()43b|) . and the uniform-electron-gas approximation for the correlation-energy 
functional, Eq. (jU^ . we obtain an approximate exchange-correlation functional, given by 

E^^'[Pi]^EM + E:M (121) 
An alternative approximation is given by the BLYP functional that uses Eqs. ()lUUp and 

dni: 

^r[Pi] ^ E:^[p] + AE^ip] + E^^p]. (122) 
If both of these approximations are reasonable, then so is a linear combination of of the two: 

^rbi] ^ b {e:^p] + AErip] + e'/j'Ip]) + (i - &) {em + em , (123) 

where the semi-empirical parameter b is at our disposal. Using the approximation for E^^, 
given by Vosko, Wilk, and Nusair jg^, denoted by we have 

^^"'[Pi] ^ b [E^^lp] + AErip] + E'y^lp]) + {l-b) {E^[p^] + E:r[p]) . (124) 

Insertion of additional empirical parameters that deviate slightly from unity, we have 

[Pi] ^ b {e:^[p] + aAE^ip] + PE'znp]) + {l~b) {EM + ^Eirip]) , (125) 
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where we can reduce the four variables to three by introducing a constraint, e.g., (7[6, /?]). 
By choosing {b = 0.80), (a = 0.90), {p = 1.01), (7 = (1 - /9&)(1 - b)-^ = 0.95) we get the 
same exchange-correlation energy functional as the one employed in the B3LYP approach. 
In particular, we have [b = 1 — Oq), {ba = a^), and {b(3 = a^), where Oq, a^, and a^, are 0.20, 
0.72, and 0.81, respectively. 
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APPENDIX A: DERIVATION OF SOME IDENTITIES 



Consider the following density functional: 

Ei[p] = [ G,{p)dr, (Al) 



and a corresponding potential, given by j^, l2l IS. 

6p{r) dp{r) 

If the integrand Gi{p) does not contain gradient terms, we can express the density functional 
as one that appears to depends on the spin-less one-particle density-matrix: 

E,[p\ = EM = j j K^i - ^2)G,{pi) dr, dr2, (A3) 

and its corresponding operator, denoted by Vi, is defined by its kernel, given by 

Operating with Vi upon an arbitrary function, 0, and using Eqs. ()A2|) and ()A4|) . we obtain 
the following: 

Vi(f){ri) = [ t;j(ri,r2)0(r2) dr2 = [ d{ri -ra) ^^ '^^^ ^ (^2) '^ra 

J J o'Pili'2,rij 

dGj (pi(ri,ri)) dGj (p(ri)) r \Ar \ 

9pi(ri,ri) dp{ri) 

Hence, the potential and operator are identical: 

Vi = v,{r). (A6) 
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Now let the energy functional have a gradient dependence: 



Ei[p\= / G',(p,Vp)rfr, 



(A7) 



where its corresponding potential is given by 



61| 



ViXr] 



6E,[p\ dG,{p) d dG,{p) 



(A8) 



Using the identity 
functional: 



5p{v) dp{r) dr dVp{r) ' 
we can express the density functional as a one-particle density-matrix 



Ei[p\ = Ei[pi] = j j 5{yi-Y2)G, (^pi(r2,ri),(Vi + V2)/9i(r2,ri)) rfridrs, 

and its corresponding operator, denoted by Vi, is defined by its kernel, given by 

5E,[pi] 



(A9) 



'^P(r2,ri) 
5{vi - rs) 



(AlO) 



dGi{pi) d dGi{pi 



d dGi{pi) 



^9pi(r2,ri) dvi 9Vipi(r2,ri) dv2 5V2Pi(r2,ri' 

Operating with -Oj upon an arbitrary function 0, using Eqs. ()A10|) and ()A8|) . and an identity, 
given by 

d dGi{p) 



d_ dGijpi) _^ _d_ dG,{pi) 

dri (9ViPi(r2,ri) dv2 dV2pi{r2,ri) J ^^^^^ dvi 9Vip(ri)' 
we obtain the following: 



(All) 



t;i(ri,r2)0(r2) dr2 

dG,{pi] 



I '5(ri - r2) ( 



dG, (p(ri)) 
dp{ri) 



d 



dGi{p^ 



d 



dGi{pi 



api(r2,ri) dri aViPi(r2,ri) dY2 5V2pi(r2,ri 



[r2) dr2 



d 



dGi{pi 



+ 



d 



dG^{pi 



dri (9Vipi(r2,ri) dr2 c)V2Pi(r2,ri 



dGdpjr,)) d dG,{p) 
dp{ri) dri dVip{ri] 



(A12) 



Hence, we have 



Vi = vdr 



(A13) 
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